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We show that, in the semiclassical hmit and whenever the elements of the Hamiltonian matrix are random enough, 
the eigenvectors of strongly chaotic time-independent systems in ordered bases can on average be exponentially 
localized across the energy shell and decay faster than exponentially outside the energy shell. Typically however, 
matrix elements are strongly correlated leading to deviations from such behavior. 
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The quantum mechanical behavior of strongly chaotic 
systems is commonly thought to be of two, apparently 
different types. On one hand, there are the time-periodic 
systems, e.g. the Kicked Rotor, for which the eigen- 
states of the one period evolution operator, U, are under 
certain conditions localized.^ When this occurs, there 
is no repulsion between most levels and the spectrum 
of quasienergies is characterized by a Poisson spacings 
distribution.'^ One the other hand, one has the time- 
independent systems, e.g. the Coupled Quartic Oscilla- 
tors, where eingenstates in a phase space representation, 
are on average homogeneously spread over the corre- 
sponding energy shell. Such states are being regarded as 
extended leading to strong level repulsion and a Wigner 
spacing distribution.^ It is the purpose of this Letter to 
show that as far as the behavior of the eigenstates and 
in particular their localization properties, the two types 
of systems are in fact quite similar and that localization 
can occur under certain conditions in time-independent 
systems as well. 

For simplicity, in what follows we shall contain our 
discussion to a subclass of the time-periodic systems, 
namely, to kicked systems of the type 



H = K{p) + V{q) J2 Ht-nT) 



(1) 



Such systems can be mapped onto Anderson models for 
the motion of electrons on a purely one-dimensional dis- 
ordered lattice. Whenever the corresponding hopping 
elements, Wn, decay faster than and the the site en- 
ergies, e„, are random enough, the resulting eigenstates 



of the Anderson model are known to be exponentially 
localized. A more direct approach is based on the fact 
that in the basis of eigenstates of K{p) the evolution 
operator, [/„,„, is a banded matrix with elements that 
have pseudo-random phases. Therefore, its eigenvectors 
are expected to behave similarly to those of a banded 
random matrix^ of the same band width, b. Since the 
latter can be analyzed using a transfer matrix formal- 
ism, the eigenvectors are exponentially localized with a 
localization length, ^, that is proportional to b^, ^ = •yb'^. 

The Kicked Rotor for example, corresponding to 
K{p) = and V{q) = kcosq, was originally intro- 

duced as a simplified but representative Poincare map of 
a time-independent system with two degrees of freedom, 
d = 2, also known as the Standard Map. It is there- 
fore natural to expect that also its quantum mechanics 
should be reminiscent of that of time-independent sys- 
tems. One obstacle that has prevented the search for this 
similarity is the fact that eigenstates of time-independent 
systems are best understood in phase space representa- 
tions through the Berry- Voros conjecture^ However, in 
order to obtain an analogous description to that leading 
to the banded evolution operator, C/„m, a basis with a 
natural ordering is required. Let, H — Hq + V, be an 
arbitrary separation of H. Using the eingenvectors of 
Hq, Vn, arranged in increasing order of the correspond- 
ing eigenvalues, Eq ^, one obtains for H a matrix rep- 
resentation, Hnm, that is banded and moreover, has di- 
agonal elements which vary on classical energy scales. 
Although the latter feature which is related to energy 
conservation is absent in the Unm matrix where diago- 
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nal elements are on average of unit absolute value, we 
find that in the semiclassical limit, h 0, the influence 
of the diagonal on the behavior of the eigenvectors is 
partially suppressed. As a consequence, in such bases, 
the semiclassical global structure of the Unm and Hnm 
matrices is very much alike. On the other hand, a sim- 
ple mechanism that generates randomness in the Unm 
matrix is absent in the case of Hnm, making the latter 
significantly less random. 

We now turn to study the structure of the Hnm matrix. 
Let us assume that all the relevant operators, namely, H, 
Hq and V, are well behaved functions of the canonical 
variables, z = (q, p) . Then, under quite general condi- 
tions, it is known^ that the energy average of the diago- 
nal elements equals to the corresponding micro canonical 
average, that is. 



where for any function, F{z), 



f dzF{z)S[Eo - Hojz)] 
f dzS[Eo - Hoiz]) 



(2) 
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Moreover, the variance of the diagonal elements, 
ajj =< H^n > — < was shown'' to be equal to 

twice the variance of the ofi^-diagonal matrix elements 
that are close to the diagonal, Cq, and correspondingly, 
cr|j = 0{h'^~^). Therefore, in the semiclassical limit the 
distribution of diagonal elements has vanishing width 
whenever d > 1. As we now show, this fact is extremely 
helpful in understanding the behavior of the off-diagonal 
matrix elements of Hnm- In the n-th row of Hnm^ the 
average distance of matrix elements from the diagonal 
measured in units of energy is 
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and using Eq. (2) 



< {AEo,nf >^ n 



2 {{Ho,H]%g} 



(5) 



where the commutator of Eq. (4) was replaced by h 
times the corresponding Poisson bracket, [...]pb. Notice 
that while Eq. (2) holds for arbitrary values of U, Eq. (5) 
only applies for small enough h. The reason is that in Eq. 
(5), the Ihs is the average of a ratio of diagonal elements 
which; in general, differs from the ratio of the averages 
on the rhs. However, if d > 1 and h 0, then the width 



of the distribution of the diagonal elements appearing in 
the denominator of Eq. (4) becomes vanishingly small, 
therefore playing the role of a constant in the averaging. 

In order to fully grasp the implication of Eq. (5), 
it is necessary to express < (Ai?o,n)^ >^^^ in terms 
of the number of states by multiplying it with the 
mean density of states, piEo), which is 0{h~'''), such 
that < (AiVo,„)2 >V2= 0{h^-'^). On one hand, for 
d > 1 and h ^ 0, < (ANo^n)'^ >^/^ diverges. On 
the other hand however, the number of states in any 
classical energy range which, in turn, corresponds to 
the size of the truncated Hamiltonian matrix. A'', is 
0{h^'^) and therefore diverges much faster. In fact, 
< (AA^o^n)^ >^/^ /N = 0{h) and accordingly, for small 
enough h, the Hnm matrix is banded. Moreover, Eq. (2) 
implies that the diagonal matrix elements, Hnn, grow on 
average as the volume of the energy shell grows. Both 
these features are absent in the traditional Random Ma- 
trix Ensembles, e.g. GOE, and this lack of structure 
leads to extended eigenvectors. A Random Matrix model 
which does include a simplified version of this structure 
is the Wigner ensemble*'^ which is composed of banded 
matrices of band width, b, with diagonal matrix elements 
that increase on average with constant rate, a. Specif- 
ically, < hnm >e— OiuSnm, whcrc < ... >e dcnotcs cn- 



semble averaging, and a 
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l + '5nm for |n — to| < b and vanishes otherwise. As in the 
case of the GOE, it is assumed that the strongly chaotic 
nature of the classical dynamics is sufficient to ensure 
that the matrix elements are uncorrelated random vari- 
ables. In what follows we shall use the Wigner ensemble 
and its properties in order to understand the behavior 
of the eigenvectors of Hnm- 

For a = the Wigner ensemble is equivalent to the 
Banded Random Matrix Ensemble (BRME) which in 
turn can be thought of as an Anderson model with ran- 
dom, long-range hopping. Similarly, it is useful to in- 
terpret the Wigner ensemble at finite a as an Anderson 
model under the influence of a constant electric field of 
strength a. This approach enables one to visualize the 
behavior of the local density of states which is closely 
related to that of the eigenvectors. In the absence of 
electric field, the Anderson model is on average transla- 
tional invariant. Therefore, the ensemble averaged local 
density of states, 



pUE, n) =< }2\{u^)n\^d{E - Ei) >e 



(6) 



where {ui)n is the n-th component of the i-th. eigenvector 

of the random matrix and Ei the corresponding eigen- 
value, is proportional to the average density of states 
itself, p{E). In particular, for the BRME, both p{E) 
and pL{E,n) are in the form of a semicircle of radius 
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2^/2b. While turning on the electric field breaks the 
translational invariancc, for small enough fields the hop- 
ping potential varies much faster than the electric one 
and the adiabatic approximation that relays on the sep- 
aration of these two energy scales is known as the slop- 
ing band picture. In this regime, one has at each site 
the same pL{E,n) as at zero field only that its center 
is shifted to include the additional electric energy that 
is increasing linearly along the lattice. This is schemati- 
cally illustrated in Fig. |l| where the allowed region inside 
the band represents the energy shell of the Wigner en- 
semble. While taking a section through the energy band 
at a fixed position, n, gives Pl{E, n), a section at a fixed 
energy is expected to give information on the behavior 
of the average eigenvector. In particular, one expects 
that the average eigenvector is constrained to lie within 
the energy shell. Independently of the sloping band pic- 
ture, the pL{E,n) for the Wigner ensemble was derived 
in Ref. || (see also Rcf. It was found that 

PL{E,n)^—f{ 7— ,g , 7 

ab ab 

where q = [a^b)^^ . For small electric fields, q ^ 1, the 
semicircle behavior persists 

f{x,q) = (4^<z)- V8g - ^2 , (8) 

and for large fields, 9 <C 1, the profile is Lorentzian, 

= -ir^, — 2 ■ (9) 

In fact, Eqs. (8-9) only hold for x <\. For a: 3> 1, / is 
the solution of an integral equation for which 

f{x,q) ~ cexp[-2a; ln(a;e~^ ^ 2q-^ ln(x/ ^))] , 

(10) 

represents an approximate solution. 

We now turn to discuss the behavior of the average 
eigenvector, that is defined as the average variance of 
the vector component at a fixed distance from the largest 
component, g{l) =< \{ui){n — nmax))\'^ >e- We find that 
the tails oi g{l) far outside the energy shell are directly 
determined by the local density of states, pL{E,n) (see 
Fig. 1^). On the other hand, inside the energy shell the 
shape of g(l) is different in the large and small q regimes. 
At large electric field, g <C 1, the disorder is too weak to 
localize the eigenvector and accordingly, like in the case 
of the tails, g{l) is determined by the local density of 
states. That is, it takes the Lorentzian form of Eq. (8) 
with I = E/a. For weak field however, g ^ 1, the band 
is only slightly sloped and as a consequence, the varia- 
tion of the local density of states is slower than the scale 



on which localization due to disorder takes place. Thus, 
the disorder is dominant and the same exponentially lo- 
calized shape as in the absence of the electric field is 
obtained (see Fig. ||). The transition between the two 
regimes is centered at qc where the hopping range is of 
the same size as the energy shell width. A rough estimate 
of qc can be obtained assuming that the width of the en- 
ergy shell does not change with a staying 4\/2& all the 
way down to qc- While such estimate gives qc ~ 0.125, 
numerically a value of qc — 0.09 is obtained. 

In order to establish the correspondence between the 
various regimes of the Wigner ensemble and the struc- 
ture of Hamiltonian matrices, one needs to use the semi- 
classical formulas for the parameters of the ensemble. 
Notice that the variances of the matrix elements in the 
Wigner ensemble are 0(1) unlike those in Hnm and 
therefore, the effective value of the electric field for the 
Hamiltonian matrix is aef — a/ao- Then, using Eqs. (2 
- 5), one obtains that q = 0{h~^) implying that in the 
semiclassical limit the Hnm matrices are in the disorder 
dominated regime and therefore exponentially localized 
eigenvectors can be found whenever the system is suffi- 
ciently disordered. This is the main result of the paper. 

Let us now study the extent to which the Wigner 
ensemble is a good model for a particular Hamiltonian 
matrix, ^2 namely, that of the Coupled Quartic Oscilla- 
tors model, 

H = + bcqt + b-^qt - a.qlql , (11) 

at be — 7r/4 and — 1.6 where the classical dynamics 
in the Poincare section appears to be fully chaotic to a 
resolution of about 0.4% of Ti which, in turn, was taken to 
be unity. Moreover, we take Hq to be composed of two 
uncoupled harmonic oscillators with frequencies wi = 
4.11 and W2 = 1.3, truncate the resulting Hnm matrix 
to the first N basis states, N — 800, and average the 
eigenvectors to obtain the corresponding g{l) function 
(see Fig. H). The comparison with the prediction of 
the Wigner ensemble that has the same average b and a 
and the same N shows large quantitative disagreement. 
In particular, the exponentially decaying shape quickly 
saturates into broad shoulders that end in a fast drop 
corresponding to the edge of the energy shell. However, 
the reasons for this discrepancy are not related to the 
intrinsic differences between evolution type matrices and 
Hamiltonian ones, namely the behavior of the diagonal 
matrix elements since these are irrelevant in the large q 
regime. Instead, it is a consequence of the fact that the 
Unm matrix of kicked systems is a lot more random, and 
thus a lot closer to the BRME, than is the Hnm matrix 
of our example to the Wigner ensemble. In the case of 
the Kicked Rotor for example, 



3 



Unm = eM-iTiTny2){-ir-^Jm-n{k) , (12) 

where J are the Bessel functions. The exponential factor 
originates from the kinetic energy term playing here the 
role of Hq and it constitutes an efficient source of ran- 
domness. In the Hnm matrix on the other hand, the Hq 
term is both additive and not exponentiated and there- 
fore, the mechanism for generating randomness in Eq. 
(12) is absent here. It is natural to expect that, in order 
to fully recover the predictions of the Wigner ensemble, 
one needs to use more complex bases than that of har- 
monic oscillators as an alternative source of randomness 
in the matrix elements of Hnm-^^ 

In summary, although eigenstates of time-independent 
Hamiltonian systems arc; indcxxl localized inside the en- 
ergy shell, the quantitative behavior is controlled by 
the degree of correlation between matrix elements. The 
study of the way in which such correlations are influ- 
enced by the choice of basis is an exciting open question 
to be addressed in future work. In contrast to the case 
of the eigenvectors, the eigenvalue statistics is greatly 
influenced by the behavior of the diagonal matrix ele- 
ments. Specifically, a finite a makes neighboring eigen- 
values correspond to eigenvectors that strongly overlap 
and the ensuing level repulsion leads to a Wigner type 
spacing distribution.^ 
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Fig. 1. The sloping band picture. The full lines corre- 
spond to the energy band edges. The dashed lines indicate 
fixed n sections which give the local density of states and 
fixed energy sections which for small q lead to the average 
eigenvector. 

Fig. 2. The numerically obtained g{l) function for the 
Wigner ensemble with a = 2 and 6 = 14 (g = 0.0179) (o) 
are compared with Eqs. (9 - 10). Moreover, the x symbols 
represent the numerical ^g{l/4:) function for the Wigner en- 
semble with a = 1 and & = 56 (same q) , verifying the scaling 
of Eq. (7) in the n direction. 

Fig. 3. The average eigenvector for the Coupled Quartic 
Oscillators (full) when the matrix has effective Qg/ = 0.013 
and b = 12.9 {q = 446) and is compared with the corre- 
sponding result from the Wigner ensemble with the same a 
and 6 = 13 (dashed). The x symbols correspond to the g(l) 
of the Wigner ensemble with the same b but a — 0, indi- 
cating that in this regime, a has almost no influence on the 
shape of the average eigenvector except for the sharp drop at 
I fs ±780 due to the band edge. 
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